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We show that within the Left-Right symmetric model, lepton number violating decays of the Higgs 
boson can be discovered at the LHC. The process is due to the mixing of the Higgs with the triplet 
that breaks parity. As a result, the Higgs can act as a gateway to the origin of heavy Majorana 
neutrino mass. To assess the LHC reach, a detailed collider study of the same-sign di-leptons plus 
jets channel is provided. This process is complementary to the existing nuclear and collider searches 
for lepton number violation and can probe the scale of parity restoration even beyond other direct 
searches. 
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The discovery of the Higgs boson [1,2] allows to test the 
mechanism of elementary particle mass generation at the 
LHC [3]. Compared to this success, the problem of neu¬ 
trino mass in the Standard Model (SM) appears acute. 
Neutrinos may be their own antiparticles [4], and lead 
to lepton number violation (LNV). The canonical way 
of searching for LNV, neutrino-less double beta decay 
(0v2/3) [5], can be induced by light Majorana neutrinos 
or by new physics [6]. The latter, needed for neutrino 
mass, can be provided by the celebrated seesaw mecha¬ 
nism [7-11]. In particular, Left-Right symmetric mod¬ 
els (LRSM) [12], designed to explain parity violation of 
weak interactions [13], embed naturally the seesaw [7, 8). 
With the left-right (LR) scale in the TeV range, 0v2/3 
may be dominated by heavy Majorana neutrino ( N ) ex¬ 
change [14, 15], which may become favored in view of the 
cosmological bound on light neutrino masses. 

A direct strategy for LNV searches at hadron colliders 
was suggested in [16] by Keung and Senjanovic (KS) [17]. 
The KS production of heavy Majorana neutrinos would 
reveal LNV and relate directly to 0u2[3 [15, 18] and lep¬ 
ton flavor violation (LFV) [19, 20]. The Dirac mass is 
predicted [21] and may be tested at the LHC through 
LNV decays, uncovering the underlying seesaw mecha¬ 
nism and connecting to electric dipole moments [21, 22]. 
Indirect constraints [23-27] played an important role and 
comprehensive analyses [28, 29] allow the LR scale well 
within the ~ 6 TeV reach of the LHC [30]. 

In this Letter we show that LHC can probe a new chan¬ 
nel, connecting Higgs physics to restoration of parity. We 
exploit the fact that the SM Higgs can have a sizeable 
mixing with the triplet that breaks spontaneously LR 
symmetry and provides a mass to N. Through this mix¬ 
ing the Higgs can decay to a pair of N [31]. Therefore, it 
can probe their Yukawa couplings via a LNV final state 
with two same or opposite sign charged leptons and four 
jets, as shown on Fig. 1. The relevant range of N masses 


* amaiezza@ific.uv.es 
t miha.nemevsek@ijs.si 
t fabrizio.nesti@irb.hr 



FIG. 1. Dominant diagram leading to LNV Higgs decay. 


typically leads to displaced vertices. Higgs decay to RH 
neutrinos was mentioned in [32] and studied in [33] with 
effective operators, pointing out the LNV character and 
vertex displacement. Here, in the LRSM, the LNV decay 
is probing the origin of N masses, just as the standard 
decays test the origin of charged fermion masses. The 
Higgs can thus act as a portal to LNV, complementary 
to 0u2/3 and the KS reaction. 

To estimate the LHC sensitivity, we implement the 
model [34], perform a simulation of the signal and the 
expected SM background, and devise cuts. To further 
enhance the search, we simulate the displaced vertices 
from N decay and highlight their importance. Given the 
current limits on the Higgs mixing [35] (see outlook [36]), 
a discovery is possible for a high LR scale, beyond the 
reach of other direct searches. 

We conclude with a discussion on alternative models 
with potential LNV Higgs decays and a short outlook on 
the related search at e + e _ colliders. 


Left-Right symmetry and Higgs mixing. Left- 
Right symmetric models [12], based on the gauge group 
SU(2)l x SU(2) R x U(1)b-l, contain a right-handed 
(RH) gauge boson W R and three RH Majorana neutri¬ 
nos. The scalar sector of the minimal LRSM [7, 8] fea¬ 
tures a complex bi-doublet $ £ (2 R ,2 R ,0) and a pair of 
triplets A l £ (3 L ,1 R ,2), A R £ (1 L ,3 R ,2): 
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In the minimal LRSM, LR symmetry is restored at 
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high energies. The potential exhibits spontaneous break¬ 
ing and since the original work [13] it has been the subject 
of several studies [25, 31, 37-39]. We focus on the mixing 
between the triplet and the SM-like Higgs and display the 
relevant terms: 


V = - - p 2 3 {A f R A R ) 

+ A (^th) 2 + p (A^ R A R ) 2 + a(d> t 'I>)(A| f A r ) . 


The trace on the parenthesis is implied and $ = 02 ^* 02 - 
The results below hold for restoration of both generalized 
parity V and charge-conjugation C [26]. 

The parameters p are fixed in the usual way, p\ = 


2Arr 


av 


Ri 


P 2 = 0, p 3 = av 2 + 2 pv R , and neutral 


scalars develop VEVs. The LR-breaking scale is set by 
( S R ) = vr and electroweak symmetry breaking is com¬ 
pleted by ( <j> 1 } = v. For clarity we stick to the case where 
(f >§ does not acquire VEV and we suppress higher v/vr 
terms. In what follows the neutral scalars <f> and S are 
the fluctuations of °) and 

Expansion of the potential (2) around the minimum 
gives the following mass matrix for </> and 5 

M 2 = 2 ( 2Av2 ° vv *) . (3) 

\avvR 2pv%) v ’ 


Its diagonalization leads to the masses of the physical 
particles to 2 ~ 4Ar; 2 — a 2 v 2 /p, m 2 A ~ 4 pv\. Here h = 
</> cos 0 — 6 sin 0 is identified with the SM Higgs boson and 
A = 5 cos 0 + 4 > sin 0 with the further neutral state. Their 
mixing angle can be large, 6 ~ (a/2 p) (v/vr). 

Since 5 is a SM singlet, mixing universally reduces the 
SM-like Higgs couplings. Recent studies [35] allow for 
sin# < 0.44 at 2er CL, nearly independently of the singlet 
mass. 


Heavy Neutrino from Higgs decay. After sponta¬ 
neous breaking, the Yukawa term £a = YnL^ArLr + 
h.c., which couples the RH leptonic doublet Lr to the 
triplet Higgs, generates the heavy neutrino mass matrix. 
This is directly proportional to the LR scale 

M n = 2 Yjv Vr , M Wr = gVR, (4) 

where g = g^.R is the SU(2)r,r gauge coupling constant. 

To probe the spontaneous origin of N mass, one should 
observe A — > NN decays and establish that Ta_> nn oc 
m 2 N . While the production of A is small due to the large 
LR scale, in the presence of 8 the gluon fusion production 
appears. Still, A may be heavy enough to be elusive, in 
fact from perturbativity of a, A < 1, one finds toa 
5TeV(0.1/sin$). More importantly, the SM Higgs can 
decay to NN. Thus the origin of neutrino masses may 
be probed by the SM Higgs boson. 

It is useful to normalize this decay rate to the leading 
SM channel h —> bb 

TjyAf ^ tan6> 2 / TOiv\ 2 / M w \ 2 f _ 4m^ 
r bb 3 \m b ) \M Wr ) V rn 2 

neglecting the bb phase space. Including the QCD run¬ 
ning and NLO corrections [40], the ratio is enhanced by 
a factor ps 2 and is shown on Fig. 2. 




FIG. 2. Decay rate of the SM Higgs h —> NN, normalized to 
the leading h —> bb channel. 


The number of N pairs produced at the 13TeV LHC 
run with 100fb _1 luminosity is simple to estimate. Tak¬ 
ing the SM Higgs gluon fusion cross-section [41] a(gg —► 
h) = 43.9 pb and Br (bb) = 57%, one gets 500 (2000) h — > 
NN events for toat = 40GeV and 0 = 10% (20%). This 
is sufficient motivation for an in-depth collider study. 

Lepton number violating Higgs decay at the LHC. 

After pair-production from Higgs decay, each N decays to 
a charged lepton and two jets via Wr, with a RH charged 
current quark flavour structure essentially identical to the 
left-handed one [42]. Due to the Majorana nature of N , 
50% of events will result in a final state of two same-sign 
leptons and four jets with no missing energy, explicitly 
signalling LNV. 

To assess the LHC sensitivity, we extend [34] the Feyn- 
Rules [43] implementation of the LRSM [44] to include 
the mixing together with Higgs gluon fusion production. 
Parton level events are simulated with MadGraph 5 [45], 
hadronized with Pythia 6 [46] and passed to Delphes 
3 [47] for detector simulation. We use MadAnalysis 5 [48] 
for cuts and event counting. Dedicated software exten¬ 
sions are implemented in each module to study the dis¬ 
placed vertices. 

The channel h —► £ ± ^ ± 4 j carries plenty of physical 
information at parton level. The total invariant mass re¬ 
constructs the Higgs mass, while the fjj invariant mass 
reconstructs the N peak. Tagging the lepton flavor iden¬ 
tifies the RH analog of the PMNS mixing matrix and 
the related Majorana mass matrix, as with the KS pro¬ 
cess [20, 49]. Notice that with such low N masses, LFV 
constraints are easily satisfied and one may expect LFV 
in Higgs decays. 

Reconstruction at detector level is more delicate. The 
Higgs is produced with a boost y(h) ~ 3 at ^fs = 13TeV 
and the N is further boosted if toat rrih/2. As a result 
the two jets from N tend to merge. In addition, the 
jets get closer to the charged lepton and typical lepton 
isolation cuts may prevent its recognition. Furthermore, 
the distribution of transverse momentum of the lepton 
(and the jets) peaks at a fairly low value of mj,/6 ~ 
20GeV. Typical detector simulation parameters forbid 
tight leptons with pr < lOGeV, causing a loss of the 
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signal by a factor of 2. Still, the N mass peak can be 
clearly observed in the /. ij invariant mass. 

The fairly long N proper decay length 


/ o \-i ^ G%m% f Mw \ 
{ n) - 16tt 3 \M Wr ) 


( 6 ) 


characteristic for this portion of parameter space [18], 
can lead to displaced N decay products. This ranges 
from sub-millimeter to a few meters, depending on mjv 
and Mw r and results in a striking LNV signature with 
two displaced vertices. 


Background estimation. Since lepton number is con¬ 
served in the SM, there is no background at parton level 
for this final state. Nevertheless, there are three distinct 
ways in which background arises: 

1. Electron charge mis-identification and secondary 
photo-production constitute a background that is hard 
to understand in absence of real data. Since at this stage 
one cannot reliably estimate this experimental effect, we 
study the muon channel free from such issues [50, 51]. 

2 . The main prompt muon background comes from 
pair-production of electroweak gauge bosons, in partic¬ 
ular WZ, ZZ and W ± W ± jj, and tt production. They 
also contain non-prompt muons from mesons. 

3. Significant background is due to non-prompt 
muons. This component, likely dominant and not easy 
to estimate, is due to QCD jets when some hadron is 
mis-identified as a muon. Although the mis-identification 
probability is small, the huge QCD production compen¬ 
sates [50, 51]. A realistic estimate will require a knowl¬ 
edge of hadron mis-id within the real detector in the next 
LHC run. Nevertheless, previous studies indicate this 
background behaves similarly to the WZ + ZZ back¬ 
ground (see supplement of [51]). From that sample, we 
estimate the QCD mis-id contribution by multiplying the 
WZ + ZZ background by 2.5. 

Selection criteria and sensitivity. We adopt the de¬ 
fault Delphes 3 ATLAS card with muon isolation param¬ 
eters in agreement with [50] and the anti-fc^ jet algorithm 
with A R = 0.4 and pTjmin = 20 GeV. We demand two 
same-sign isolated muons and no other leptons, together 
with n.j jets, where rij = 1, 2,3. We require $ T < 30 GeV 
and leading muon transverse momentum px < 55 GeV. 
We demand the transverse mass < 30 GeV and in- 

variant masses m < 80 GeV, < 60 GeV. The 

impact of these selection cuts is shown in Tab. I. 

For both short and long lived Ns the known decay 
length allows us to impose cuts on the muon vertex trans¬ 
verse displacement dx, shown on Fig. 3. We smear the 
reconstructed vertex with the px—p dependent resolution 
of 20-40 pm, as reported in [52]. Since the typical back¬ 
ground contains one prompt and one secondary muon, it 
is effective to cut on both short and long dx, on both 
muons. We employ a sliding window cut, allowing events 
with L /10 < dx < 5L and optimizing over L. This cut 
is expected to give further control on the multijet QCD 
background. 


Process 

No cuts 

Imposed cuts 

p^+rij 

J5 T 

Pt 

nriT 

777-inv 

WZ 

2 M 

544 

143 

78 

40 

20 

ZZ 

1 M 

55 

29 

16 

12 

8 

W ± W ± 2j 

389 f 

115 

16 

5 

3 

1 

it 

10 M f 

509 

97 

40 

22 

14 

Signal (40) 

543 

44 

43 

41 

38 

37 


TABLE I. Number of expected events at the 13TeV LHC 
run with C = 100fb _1 after cuts described in the text. The 
signal is generated with 40 GeV, sin 6 = 10%, Mw R = 3TeV 
and rij = 1, 2, 3. ^Here we restrict to t —> bW and W —» pu. 


The final sensitivity is shown in Fig. 4, where a single 
N with a 100% Br to muons is assumed. Given the al¬ 
lowed amount of Higgs mixing sin 9 < 40%, this channel 
allows to test the LR scale beyond the expected reach 
Mw H ~ 6TeV of direct searches [30]. 

Other models with LNV Higgs decays. The heavy 
fermions S in standard type-I [7-11] or type-III [58] see¬ 
saw couple to the Higgs boson via their Dirac masses 
rriD- If these were rather large (unlikely in LRSM [21]) 
h —► vS decay might be observed [59] and via further 
mixing mo /ms, the LNV channel h — > SS could open 
up [32]. Below the Z mass, where such two body Higgs 
decay is available, constraint on the mixing [51, 60] sup¬ 
presses it beyond observation. 

Further extension with a singlet [61] or spontaneously 
broken B — L models can generate an observable signal. 

Models with supersymmetric R-parity violation, alter¬ 
native to the seesaw mechanism [62], lead generically to 
LNV. Same-sign di-lepton decay of the Higgs boson could 
follow from a fairly large h —> x°X° coupling, where 
decays to ijj thanks to a non-zero Iqq' term. This sce¬ 
nario may deserve an updated study in light of recent 
limits. In case of slepton-sneutrino mass degeneracy, the 
LNV mode appears to be suppressed [63]. 



FIG. 3. Reconstructed transverse muon displacement after 
g ± g ± +rij event selection and before other cuts. 
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FIG. 4. Expected sensitivity in the Mw R — rriN plane, in the 
cases of small and large allowed mixing sin (9. The current 
bound from the KS channel [53] is shown in the upper left 
corner; the limit from di-jet search [54] as the vertical band. 
The region excluded by the GERDA I [55] 0u2/3 search is the 
diagonal shaded region with a factor of 4 uncertainty from 
matrix elements [56], while the red line marks the sensitivity 
of GERDA II. The grey dashed lines mark the lab frame decay 
length of N from Higgs decay and the lower golden triangle 
the region where the length from Wr —> IN decay exceeds 
5 m, enabling W' — > 1$ searches [57]. 


dent measurements. The event rate and invariant mass 
peaks give information on 6 x Y}v and tun , respectively. 
The decay length then constrains AT Wr via Eq. (6) and 
finally a global fit on the Higgs data can provide sin#. 
If more than one N were observed, this task simplifies. 
To complete the understanding of neutrino mass origin, 
one would clearly like to observe A and the associated 
gauge boson that provides the gauge symmetry protec¬ 
tion. This channel would nevertheless allow to probe the 
polarization of N decay, similarly to the KS channel [64], 
even if Wr were out of reach. A comprehensive study 
will be presented in [65]. 

A number of potential improvements can be identified: 
i) adding muons with px < lOGeV would double the 
signal; ii) for short lived N, tight cuts on displacement 
could reduce further the QCD multijet background; iii) 
for boosted long lived Ns, the muon tends to merge with 
the jet and one gets displaced jets. Although a challenge, 
it may be possible [66] to identify their displacement. 

We conclude by pointing out that e + e _ colliders pro¬ 
vide a particularly clean environment for heavy neutrino 
searches [67]. For the LNV Higgs decay, the relative 
decrease in production cross-section to cr = 0.24 fb at 
yfs ~ 240 GeV may be compensated by lack of back¬ 
ground (only ZZ remains) and large luminosity ~ 1 — 
10ab” 1 2 3 [68]. Conversely, a positive signal at LHC with¬ 
out the associated Wr discovery would make a case at a 
high energy hadron collider. 
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